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Differential stability of DNA based on salt concentration 
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Intracellular positive ions neutralise negative charges on the phosphates of a DNA strand con¬ 
ferring greater strength on the hydrogen bonds that connect complementary strands into a double 
helix and so confer enhanced stability. Beyond a certain value of salt concentration, the DNA 
molecule displays a unstable nature in vivo as well as in vitro. We consider a wide range of salt 
concentrations and study the stability of the DNA double helix using a statistical model. Through 
numerical calculations we attempt to explain the different behaviour exhibited by DNA molecules 
in this range. We compare our results with experimental data and find a close agreement. 


I. INTRODUCTION 

Intracellular sodium triggers a cell to progress to cell 
division [^. In many biological activities like creation 
of daughter cells, salt enacts a pivotal role. It is also a 
known fact that the two strands of DNA carry negative 
charge due to their phosphate groups {POjf ) 0. The two 
strands are connected with each other through the hy¬ 
drogen bonds between the bases on the opposite strands. 
However, the stability of the DNA double helix also de¬ 
pends on the charge neutrality of the system. To counter 
the negative charges on the strands, the cellular environ¬ 
ment requires some positive charges in order to nullify 
the repulsion between negative charges. Such positive 
counter ions derive from salts present in the body. The 
primary salts are NaCl and MgCl 2 which on reaction 
with water release iVa+ and while the chlorine 

rearranges itself in the form of hydrogen chloride (HCl). 
These cations act as electrostatic screening agents by cre¬ 
ating an attractive force between the positive ions [1, Q 
and the negative charges on the strands. The stability of 
the DNA helix is then governed by the balance between 
the attractive and repulsive forces exist in the system. 
These forces exist in the form of hydrogen bonding be¬ 
tween the bases on the opposite strands, the stacking in¬ 
teraction along the strands, the repulsion between same 
kind of charges and attraction between opposite kind of 
charges. Various experiments followed by theoretical in¬ 
vestigations on the role of salt in screening the repulsive 
forces in DNA revealed many interestiM results. Not 
only the stability of DNA molecule pl-Q but effects of 
salt on the B — A transition, and on the condensation of 
the DNA molecule have been studied in detail by several 
researchers In all these studies the stability of DNA 
has been studied under low or moderate (O.I-I.O M) con¬ 
centrations of salt. Several exp eriments, conducted by 
the groups of Owczarzy [13, [Hf and SantaLucia [l^, [HI 
revealed that short as well long DNA molecules become 
more stable as the content of cations in the solution in¬ 
creases, with the stability of the structure being deter¬ 
mined by calculating the melting temperature (T^) of 
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the molecule. In all these experiments, the melting tem¬ 
perature was found to have a logarithmic dependence on 
the amount of cations present in the solution. The the¬ 
oretical work following these experiments attempted to 
explain this logarithmic dependence, based either on a 
semi- emp irical formula [3, [iJ, ll3 or on statistical models 
d, ITbHl^ . The counterion theory proposed by Man¬ 
ning explains the stability of the DNA molecule due to 
the presence of salt in the solution Bi- The stretch¬ 
ing and unzipping behaviour of the DNA molecule in the 
presence of cations have also been discussed by several 
groups B [l^ 113 . These results show that the mechani¬ 
cal stability of the DNA molecule also increases with the 
concentration of cations. 

In another set of experiments [llj, [13 that were exe¬ 
cuted at relatively high salt concentration, some strik¬ 
ingly different behaviour in the DNA molecule was ob¬ 
served. These experiments found that in this range of 
salt concentrations the stability of DNA gets shattered. 
Interestingly, in a similar set of experiments, several re¬ 
search groups tested the condensation process and melt- 
in g pr ofile in DNA under different ethanol concentrations 
[sL l23| . One of the findings of these experiments revealed 
that at a relatively high concentration of ethanol, the 
melting temperature of the system increases. The role 
of cations in vivo is manifold. In a recent review, Bose 
et al [13 showed that abundant cations inside the cell 
as well as outside the cell play important roles in sus¬ 
taining cancer cells and at the same time in the decay of 
immune cells. Cancer cells possess a significant electrical 
character [13 compared to normal somatic cells. Cancer 
cells become more electrostatic [l| during oncogenesis. 

Keeping in view the important aspect of salt in the cell, 
we propose a theoretical description of the stability of 
DNA molecules at higher salt concentration. We consider 
one of the DNA sequences studied by Khimji et al. [13, 
who studied the DNA duplex stability in crowded polyan¬ 
ion solutions. For the current investigations we adopt the 
standard Peyrard Bishop Dauxois (PBD) model [H, [13 
and modify the potentials appearing in the model. In 
our earlier work, we showed that the PBD model has 
enough detail to explain the thermal as well as mechani¬ 
cal stability of DNA molecules at lower strengths of salt 
and with varying salt concentration B, [^ ■ Our presen¬ 
tation in the current manuscript is organized as follows. 
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In section im we discuss the PBD model and the suitable 
modifications we did to incorporate the high concentra¬ 
tion of cations. How the stability of DNA in thermal as 
well as in force ensembles varies in a wide range of salt 
concentrations in solution is discussed in section mi We 
finally conclude our results in section ITVl 


where 7 is the solvent interaction factor and it reduces 
the height of the barrier appearing in the potential 
The concentration of salt that stabilizes the DNA 
molecule can be incorporated through the on-site poten¬ 
tial term in the Hamiltonian. The presence of 

cations around the DNA molecule are shown, schemat¬ 
ically, in Fig. [H It is a known fact that the cations 


II. MODEL 


We adopt the PBD model to study the response of 
DNA in different ensembles. We use the linear version of 
the PBD model. Although this version simplifies the he- 
licoidal geometry of the molecule, it has enough details to 
delineate the stability of the molecule in thermal as well 
as in force ensembles [ 2 ^-[ 3 l|. The model represents the 
hydrogen bonding between the bases in a pair through 
the Morse potential while the stacking interaction is rep¬ 
resented by an anharmonic potential. The Hamiltonian 
of the system is. 
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^ = E 
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N-1 

+ E 

2=1 


( 1 ) 

where {pi = myi) represents the momentum of a base 
pair where m is the mass of a base pair and yi represents 
the stretching of the base pair. The stacking inter¬ 
action between the nearest neighbours is represented by 
an anharmonic term, Vs{yi,yi+i)^ which is. 


k 

Vsiy^,y^+l) = -{yi - y^+i^ll + pexp{-b{y, + yi+i)}]. 

( 2 ) 

Here k represents the single strand elasticity. The an- 
harmonicity in the strand elasticity is taken care of by 
the term p while b represents the range of anharmonic- 
ity. For our studies we tuned the model parameters and 
found the values of fc = 0.01 ev/A^, p = 1.0, b = 0.035 
A as suitable values. The potential Vm{yi) describes the 
interaction between the two bases in a pair [2^, 

VM = - 1 )" ( 3 ) 


where Di is the equivalence of the dissociation energy of a 
pair. The constant, a^, represents the inverse of the width 
of the potential well. These two parameters have a crucial 
role in DNA denaturation. The dissociation energy is a 
representation of the hydrogen bond energy that binds 
the A — T and G — C. It is a known fact that the bond 
strength of these two pairs are not the same but are in an 
approximate ratio of 1 . 25 - 1.5 as the GC pairs have three 
while AT pairs have two hydrogen bonds [Tsl. I 3 l - l 4 ll| . The 
last term in the Hamiltonian is the solvent term Vsoi{y)- 
This is an additional term which simulates the formation 
of H-bonds with the solvent, once the hydrogen bonds are 
stretched by more than their equilibrium values. This is 
expressed as[l, 

Vsoiiyi) = -jD.ltanh^-fy,) - 1 ] ( 4 ) 



FIG. 1. The schematic representation of cations and DNA. 
The figure shows the probable distribution of cations around 
the two negative strands which are standing around 20A apart 
and the cations trying to nullify the repulsion between the 
two strands and possible distribution of excess cations that 
interact among themselves. 


are required to reduce the repulsion between and along 
the DNA strands so that the DNA is stabilised in its 
double stranded configuration. However, there is a criti¬ 
cal concentration of cations in the solution below which 
these cations act as shielding agents. In case the num¬ 
ber of cations in the solution crosses this critical num¬ 
ber, the repulsive forces in the system dominate over 
the attractive forces which further shatter the stability 
of the dsDNA. The overall activities in the surround¬ 
ings of DNA molecule at higher concentration might be 
complicated to analyse. However, as an intuitive guess, 
one can anticipate the behaviour of cations as Coloumbic 
particles. We are helped by the basics of electrostatics 
here. From elementary knowledge we know that for a 
monopole the potential scales as i, where r is the dis¬ 
tance between the charges. Similarly for dipole it scales 
as for quadrupole ^ and for octupole as ^ etc [i^ . 
We extend this concept and introduce a scaling term that 
represents the repulsive contribution in the complete po¬ 
tential. The cation’s concentration, C around the DNA 
molecule can be mapped with the inter particle distance 
r. The logic behind this proposal lies in the fact that 
when C increases, r will decrease and vice versa. In the 
PBD model, the prime contributor to the stability of the 
molecule is the potential depth of the Morse potential 
{Di). Thus, we adopt this parameter {Di) as a function 
of the salt concentration of the solution. Considering all 
of the above forces that are responsible for the stability 
of DNA molecules, we express the parameter, Di, as (in 
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terms of salt concentration), 


Di — Do 


1 + Ai In 






,(5) 

where the first three terms are taken from our previous 
work 0 which show the logarithmic dependence of Tm 
on salt concentration. Here we make (7 as a dimension¬ 
less quantity by choosing the reference value Cq = 1.0 
0 . The fourth term is added which takes care of the be¬ 
haviour of cations in the surrounding of DNA molecules 
at high salt concentrations. This is a pure Coloumbic 
term. A point to note is that t is an important param¬ 
eter. Since the nature of cations and their exact dis¬ 
tribution is not known, how the Coloumbic interactions 
modify the interaction between the two strands is not 
precisely known. At the same time, it is also a fact that 
cations shield the repulsion between the negative ions. 
We are not sure about the nature of the pole that might 
be existing due to presence of these cations. The best 
we can think of is that they obey some power law. We 
propose that the potential depth will scale with the con¬ 
centration of salt as C~*. This term will have a pro¬ 
portionality factor X which can be found by comparing 
the theoretical results with the experimental findings. At 
lower values of C, the logarithmic term will dominate, so 
the molecule will be more stable with increasing values of 
C. For higher concentration the interaction between the 
cations will dominate. For our calculations we found the 
best match for t = 0.01 and y = 1.2. The other potential 
parameters are Dq = 0.043 eV, Ai = 0.01, A 2 = —0.011, 
OAT = 4.2 A acc = 6.3 A We take 7 = 1.0 A ^ in 
this problem as a suitable choice for our calculations. 

The partition function can be written with the com¬ 
plete Hamiltonian of the system like this, 


Z = 


N 

Y\_exp[-PH{yi,yi+i,pi)]dyidpi ( 6 ) 


where /3 = -j^^- The momentum part can be easily in¬ 
tegrated and is equal to (27rmfcBT)^/^. The configura¬ 
tional partition function is defined as, 


N 

= / Y[K{y^,yi+i)dy^ (7) 

i—1 

For the homogeneous chain, one can evaluate the parti¬ 
tion function by transfer integral (TI) method by apply¬ 
ing the periodic boundary condition |44| . In the case of a 
heterogeneous chain, with an open boundary, the onfigu- 
rational part of the partition function can be integrated 
numerically with the help of the matrix multiplication 
method (TII - It^ . The important part of this integration is 
the selection of proper cut-offs for the integral appearing 
in Eq .6 to avoid the divergence of the partition function. 
The method to identify the prop er cut-off has been dis¬ 
cussed by several groups [13, Es Eli ■ In the calculations 


by T.S. van Erp et al it was shown that the upper cut-off 
should be ~ 50 A with our model parameters at T = 300 
K while the lower cut-off is -0.3 A. In the earlier work by 
Dauxois and Peyrard it was shown that the Tm converges 
rapidly with the upper limit of integration (^ . In this 
work the equation for the partition function for an infi¬ 
nite homogeneous chain was solved using the TI method. 
For a short chain, we calculated Tm for different values 
of upper cut-offs and found that an upper limit of 200 A 
is sufficient to avoid the divergence of the partition func¬ 
tion. Thus the configurational space for our calculations 
extends from -5 A to 200 A. Once the limit of integra¬ 
tion has been chosen, the task is reduced to discretize 
the space to evaluate the integral numerically. The space 
is discretized using the Gaussian quadrature formula. In 
our previous studies (4^ . we observed that in order to 
get precise value for melting temperature (Tm) one has 
to choose large grid points. We found that 900 is a suf¬ 
ficient number for this purpose. As all matrices in Eq. 
[7] are identical in nature the multiplication is done very 
efficiently. The thermodynamic quantities of interest can 
be calculated by evaluating the Helmholtz free energy of 
the system. The free energy per base pair is, 


fiT) 
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( 8 ) 

The thermodynamic quantities like specific heat {Cy) as 
a function of temperature or the applied force can be 
evaluated by taking the second derivative of the free en¬ 
ergy. The peak in the specific heat corresponds to the 
melting temperature or the critical force of the system. 

Other quantities like the average fraction 9{= 1 — </>) of 
bonded (or open) base pairs can be calculated by intro¬ 
ducing the dsDNA ensemble(dsDNA) or using the 
phenomenological approach [4^, [^. In general, the 0 is 
defined as. 


0 = 0e.t0int (9) 

0 ext is the average fraction of strands forming duplexes, 
while 0 int is the average fraction of unbroken bonds in 
the duplexes. The opening of long and short chains are 
completely different. For long chains, when the fraction 
of open base pairs, ^(1 — 0), goes practically from 0 to 
1 at the melting transition, the two strands are not yet 
completely separated. At this point, a majority of the 
bonds are in a broken state and the dsDNA is denatured. 
However, few bonds are still in an intact state, prevent¬ 
ing the two strands separating from each other. Only at 
high temperatures there is a real separation. Therefore 
for very long chains the double strand is always a sin¬ 
gle macromolecule through the transition, thus one can 
calculate the fraction of intact or broken base pairs only. 
While for short chains, the process of single bond disrup¬ 
tion and strand dissociation tends to happen in the same 
temperature range. Thus, the computation of both 0mt 
and 0ext is essential [50j |. 






III. RESULTS 


340 


First we discuss the thermal denaturation of the DNA 
molecule. To find the thermal stability of the molecule 
we calculated the melting temperature with the help of 
Eq. IHiyni Here we discus the stability of the short DNA 
sequences that are used by Khimji et al [l^. We consider 
the sequence 5' — TCACAGATGCGT — 3' for our stud¬ 
ies. The peak in the specific heat and in the differential 
melting curve indicate the melting temperature of the 
molecule (see Fig. [2]). We calculate the melting temper¬ 
ature ~ 

The ] n n 



FIG. 2. The specific heat as a function of temperature 
for the 12 base pair chain for the model parameters k = 
0.01 eV/A^, Do = 0.043 eV, Ai = 0.01, A 2 = -0.011, t = 
0.01, X — 1-2. The figure is plotted for C = 0.5 M. The peak 
corresponds to the melting temperature. 



FIG. 3. The temperature-salt phase diagrams showing the 
variation of Tm with salt concentration. The comparison of 

01 n 



the figure we find that the first three terms in Eq. (4) are 
responsible for the behaviour of the molecule in the low 
salt region, between 0.01-1.0. Once the concentration is 
above 1.0 the Coloumbic repulsion among the cations in¬ 
creases. This leads to the instability in the molecule. As 
discussed earlier, the parameter t plays an important role 
here. We calculate the phase diagram for different values 
of t and find the best match for t = 0.01 for the exper¬ 
imental results of Khimji 22|. The maximum deviation 
is at G = 5.0 which is about 0.8 K. Thus, we can say 
that our results are in close match with the experimental 
findings of [^ . 

In order to validate our results for long chains, we in¬ 
creased the length of the molecule by repeating the same 
sequence. From Fig. [4] we find the behaviour of phase 
diagram is the same with increasing lengths of DNA 
molecule. However, as the chain length increases the 
phase diagram for different lengths comes closer. This 
indicates the limit of infinite chain length. We do not 
have the experimental data for infinite chains at this mo¬ 
ment to compare our results. 

It is important to note here the difference in the na¬ 
ture of the slope at lower and higher concentrations. Ex¬ 


FIG. 4. The phase diagram for different chain lengths in order 
to validate the results for infinite chain. The lengths are 36, 
60, 156, 252, 276 bps. 


perimental as well as theoretical findings clearly indicate 
that the stability rate is much sharper than the instabil¬ 
ity rate. The prime reason behind this difference might 
be the effect of excluded volume at higher concentrations. 
As the concentration of the salt increases the number of 
cations increase correspondingly and therefore the vol¬ 
ume available for the movement will lessen. The cations 
have charge as well as size and so occupy the flexible 
volume (they crowd the DNA). At higher concentrations 
although there is a instability due to repulsive forces, 
there is a stability due to excluded volume. While the 
movement of cations at lower concentration is faster, at 
high concentration their movement might be hampered. 

After evaluating the stability of DNA in a thermal 
ensemble we calculated the stability of the molecule in 
a force ensemble for the above mentioned range of salt 
concentration in the solution. It is a known fact that 
the replication process is initiated by the force exerted 
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by DNA polymerase on a segment of DNA chain. Repli¬ 
cation starts at the replication origin [Ell, and the 
replication fork propagates bidirectionally. Mathemati¬ 
cally one can model the replication as the force applied 
on an end of the DNA chain [S^IU. The physics of open- 
ing the chain in these ensembles is completely different 
(2Cl ISSjl .Whereas in thermal denaturation, the opening is 
due to an increase in the entropy of the system, for a me¬ 
chanically stretched DNA chain the opening is enthalpic. 
The modified Hamiltonian for the DNA that is pulled 
mechanically from an end is, 


N r 






2m 




N-1 


E {^s{yi,yz+i)]-F 


i=l 



thermodynamic quantities of interest, from the modified 
Hamiltonian, can be calculated using Eq. 0 & [5] Here 
we consider an infinite chain of 300 base pairs which is 
generated by repetition of the 12 base pair chain [2^ . 
In vitro, experiments on DNA unzipping are executed 
either at constant displacement or at varying loading 
rates. Although in both the set-ups, the microscopic and 
dynamic behavior of unzipping are different, the criti¬ 
cal force comes out to be the same for an infinite chain 
(or in thermodynamic limit), i.e. for A-phase DNA. For 
theoretical investigation on force induced unzipping of 
DNA, the chain length of the molecule should be large 
enough to consider the chain as infinite. We calculated 
the melting temperature, Tm for different chain lengths 
(in increasing order) and found that a length of about 
300 base pairs is sufficient to consider the chain as an 
infinite chain [E^ . All the base pairs of dsDNA kept in a 
thermal bath share equal amount of energy while in the 
case when the molecule is pulled from an end there is 
a differential distribution of the applied force (from the 
pulling point to the other end of the chain). The phase 
diagram for the mechanically pulled DNA shows that the 
force required to open the chain decreases at higher salt 
concentrations in the solution, as shown, FigjSl Since 
the molecule itself is unstable at higher salt concentra¬ 
tions, the force required is less. Unfortunately, we do not 
have any experimental data for force induced unzipping 
at higher salt values. We hope that some experiments 
may be performed to validate these results. 



FIG. 5. The critical force-salt phase diagram showing the 
variation of iT with salt concentration. As there is no ex¬ 
perimental data is available to compare, we show only the 
calculated points on the graph. 


cations act as shielding agents for the negatively charged 
strands of DNA. However, there is a range up to which 
they contribute to the stability of the molecule. When 
the concentration exceeds this range, their role in sta¬ 
bility may be shattered [E^. Our work is an attempt 
to explain the possible mechanism that might be taking 
place at higher concentrations of salt in the solution. The 
volume available in the cell is fixed and there is a limit 
up to which these cations contribute to the overall charge 
neutrality of the system. A greater number of cations in 
a fixed volume may disrupt the balance between the pos¬ 
itive and negative charge forces. Through the schematic 
diagram and by calculating the free energy as a function 
of temperature and force, we explain the cause of the in¬ 
stability of the DNA molecule at higher concentrations. 
Our results are in very good agreement with the exper¬ 
imental findings [E^. However, due to unavailability of 
results for infinite chains we can not compare our results 
for force induced unzipping of DNA molecules at high salt 
concentration. It would be interesting to calculate the ki¬ 
netics of duplex unzipping at higher salt concentration. 
The movement of cations at higher salt concentrations 
might also be interesting to analyse. 


IV. CONCLUSION 

In the current work, we propose a theoretical descrip¬ 
tion of the behaviour of DNA molecules at relatively 
high concentrations of salt. There are several biologi¬ 
cal events where the concentration of cations increases in 
the cell. For example, in somatic cells the concentrations 
of sodium ions substant ially increase during activation to 
become cancer cells [I|, [E2|- From most of the previous 
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